Effect of optical phonon scattering on the performance of GaN transistors 
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A model based on optical phonon scattering is developed to explain peculiarities in the current 
drive, transconductance, and high speed behavior of short gate length GaN transistors. The model 
is able to resolve these peculiarities, and provides a simple way to explain transistor behavior in any 
semiconductor material system in which electron-optical phonon scattering is strong. 
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GaN-based high electron mobility transistors 
(HEMTs) exhibit a number of peculiar features. Polar- 
ization induces some of the highest channel 2D electron 
gas (2DEG) sheet charge densities in all semiconductors, 
with values in the range of n s ~ 1 — 4 x 10 13 /cm 2 . If an 
effective electron saturation velocity of v sat ~ 10 7 cm/s 
is assumed, the drain current per unit device width 
(W) is expected to be I s d at /W = J sat = qn s v sat ~ 5 
A/mm for n s ~ 3 x 10 13 /cm 2 . However, experimentally 
measured saturation currents fall far short of such values 
- which is the first peculiar feature. Second, the 'gain' 
of GaN HEMTs, or their transconductance given by 
g m = dId/dVg S shows a sharp drop after reaching a 
maximum as the gate voltage is increased. This is also 
not expected, since the constant saturation velocity 
model predicts g m /W = Cg S v satl which should stay 
constant at high gate voltages. Third, the current-gain 
cutoff frequency (fr) also follows the trend of g m . 
Earlier explanations attribute this to a variable source 
resistance [1]. However, recent results [2] with negligible 
source access resistances (~ 50 Q - fim) show an even 
pronounced 'peaky' behavior of g m and fr vs V gs . The 
purpose of this work is to offer a novel, yet simple model 
of GaN transistor operation that reconciles all three 
peculiarities. In addition, the model also brings forth a 
remarkably simple picture of transistor operation when 
electron-optical phonon scattering is strong, one that can 
be used for quantitative predictions. It also highlights 
the importance of an accurate representation of the 
current drive in a transistor, a central quantity that 
determines all other device performance characteristics. 

The first hint at resolving the above peculiarities 
appears upon a close examination of the strength of 
electron-optical phonon interactions in GaN. Owing to 
the light mass of the nitrogen atom (lightest in all III- 
V semiconductors), the polar optical phonon energy is 
high. Combined with the high electronegativity of nitro- 
gen, it results in a very high scattering rate. The mean 
free path of energetic electrons emitting optical phonons 
is X op ~ a^e oc /(eo — e^) ~ 3.5 nm (a^ is the effective 
Bohr-radius and eo, are the static and high-frequency 
dielectric constants). This is much shorter than most 
other semiconductors (Si, GaAs). It prevents ballistic 



transport at high- fields, and damps velocity overshoot 
effects even for very short gate length GaN devices. Sim- 
ilar hi-field phenomena also occurs for other high phonon 
energy materials such as carbon nanotubes and graphene 
[3-5]. 

Fig 1 depicts the band diagram of a typical GaN tran- 
sistor. To understand the characteristics of a short gate- 
length device, we inspect the electron distribution at 
the 'source injection point' [6]. The energy and the 
k— space distribution of 2DEG electrons at this point 
is depicted for various carrier densities dictated by the 
gate voltage and the gate-to-source injection plane ca- 
pacitance. We assume an electrostatically well-designed 
device with negligible drain-induced barrier lowering and 
other short channel effects. Fig l(a, b, and c) depict 
the carrier distributions at equilibrium (Fermi circles of 
radius kp = \Z2i\n s at Vds = 0, n s is the 2d sheet den- 
sity), and at current saturation (shifted Fermi circle of 
the same radius). In a ballistic transistor [6], the differ- 
ence of quasi-Fermi levels of the right- and the left-going 
carriers is the drain bias. However, in GaN transistors, 
ultrafast optical phonon emission locks this difference at 
the optical phonon energy (Hlj op ~ 92 meV), as shown in 
Fig 1. The carrier distribution in k— space at current sat- 
uration is depicted by the filled circles at various 2DEG 
carrier concentrations; as the gate pinches the channel 
off at the source injection point, the Fermi circle shrinks. 

At high 2DEG densities (Fig 1(a)), electrons from the 
highest right-going energy state emit an optical phonon 
and scatter into the highest empty left-going state. Scat- 
tering into the bottom of the subband is Pauli-blocked 
due to the high degeneracy, but becomes possible when 
the condition 2kp < k op is met as shown in Fig 1(b). 
From Fig 1(b), we get H 2 kl p /2m* = Huj op , and the 
crossover 2DEG density is n = k^ /Sir ~ 1.9 x 10 12 /cm 2 
(here m* ~ 0.2mo is the electron effective mass of GaN). 
For n s < no, the highest energy right-going state has an 
energy hw opi and the lowest occupied state also moves 
to the right, as seen in Fig 1(c). Now that the carrier 
distribution at the source injection point in the k— space 
is identified, the current flowing in the HEMT can be 
calculated. The current is found by summing over the 
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FIG. 1: Transistor cross section, energy band diagram, band- 
structure, and k— space occupation of electrons. 

group velocities of all occupied states in the k— space; for 
2D carriers in a parabolic band, the result is simple: 

J =7^W **(^)/00 = (1) 
{Z7r) z J r m* m* 

where k is the geometrical centroid of the occupied area 
in the k— space. The result is general, valid at all tem- 
peratures if /(k) is symmetric around ko, and is not re- 
stricted to GaN transistors. For the specific case con- 
sidered here, k = kox, and the scalar drain current is 
J = qn s (hko/m*). Thus, the net current per unit width 
under any bias conditions can be found if the centroid of 
the distribution in k— space is identified. 

For high sheet densities (n s > no ~ 1.9 x 10 12 /cm 2 as 
depicted in Fig 1(a)), the centroid is found by solving 

h 2 (k F + k ) 2 _ h 2 (k F -k ) 2 = 

2m* 2m* op ' U 

yielding fco = m*uj p/2hkF. Similarly, for the low-carrier 
density condition in Fig 1(c), it is easily seen that the 
centroid is at fco = k op — hp. The HEMT saturation 
current density is thus given by 



J= qn a n ^-J 2 ™'> n s <n , (3) 
= ^cj opv /n~ n s > n . (4) 

The results indicate a significant deviation from a con- 
stant saturation velocity model. The current is propor- 
tional to the square root of the 2DEG density in the high 
density (and high current) regime. The difference is high- 
lighted in Fig 2(c). 



The transconductance in S/m units is given by g m = 
dJ/dVg = (dJ/dn s )C gsi where C gs is the intrinsic ca- 
pacitance per unit area between the gate metal and the 
source injection point. Using Eqs. (3, 4), we get 

gm = ^r(ko P ~ ^V 2 ^s)C gs n s <n (5) 

= i^fe:^ > "o> ( 6 ) 

which casts the transconductance in a familiar form as 
the product of the gate capacitance and the prefactor 
indicative of an effective electron ensemble velocity. Note 
that the effective velocity is strongly dependent on the 
carrier density, and therefore we are forced to refrain from 
invoking the concept of a saturation velocity that is a 
material constant. Furthermore, it is evident that the 
transconductance at high carrier densities will decrease 
as g m ~ \j yfn~ s . With a parasitic source resistance R s 
depicted in Fig 1, there is a voltage drop from the source 
contact to the source injection point, and the measured 
transconductance becomes g^ 1 = g m /(l + R s gm)- 

The speed of the transistor can be evaluated by either 
linearizing the small-signal temporal drain current (Eqs. 
3, 4) response to a sinusoidal gate voltage sweep, or by 
a straightforward charge control analysis; both methods 
lead to the same intrinsic channel delay. The delay is 
r c = dQ/dJ = qLgdn s /dJ, which yields 

= 4 ~^F L 9 »- > n . (8) 

We note that the intrinsic channel delay is of the form 
r c = Lg/v e ff, where the effective velocity is strongly de- 
pendent on the 2DEG carrier concentration controlled by 
the gate voltage. The unity current-gain cutoff frequency 
of the transistor is then calculated as fa = l/27TT tot , 
where r tot = T c +C gd / g m +Cg d (R s +R d ), where C gd is the 
gate-drain overlap capacitance, and R s , R d are the source 
and drain resistances (units- ^.mm). All these lumped 
elements are depicted in Fig 1. The transconductance is 
the 'intrinsic' g m from Eqs. 5, 6. The dependence on g m 
indicates that in the high density regime, fa will decrease 
with increasing channel charge or equivalently, the gate 
voltage. 

The transistor characteristics under current satura- 
tion are now captured by Eqs. 3-8. The only pa- 
rameter necessary to quantitatively investigate particu- 
lar cases is the gate-channel capacitance C gs . To model 
a realistic GaN HEMT, we use a self-consistent Poisson- 
Schrodinger solver to calculate C gs as a function of V gs . 
The calculation gives an accurate dependence of charge 
at the source injection point on the gate voltage. The 
calculation accounts for quantization effects in the chan- 
nel 2DEG. The device structure chosen for highlighting 
the phonon saturation model is an Alo.35Gao.65N/GaN 
HEMT with a £& ar = 5 nm barrier. Fig 2 shows the 
device characteristics using Eqs. 3-8 using the self- 
consistent Poisson-Schrodinger results as the input. The 
calculated C gs is shown as a function of the gate voltage 
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FIG. 2: AlGaN/GaN transistor characteristics calculated us- 
ing the phonon model. 

in Fig 2(a), indicating a threshold voltage of ~ —1.5 V. 
This figure also shows the 2DEG density at the source 
injection point for three different R s values. This sheet 
density depends on the potential difference between the 
gate and the source-injection point (V g ° s ). For R s ^ 0, 
the external gate voltage is V gs = Vg S + JR S1 and this 
is what is plotted in Fig 2. As the source-resistance in- 
creases, a higher gate voltage is necessary to induce the 
same 2DEG charge. The specific case for n s = no is high- 
lighted by a red line and three arrows in Fig 2(a). This is 



also the 'crossover' density at the source injection point 
as identified earlier. The 2DEG density at the source in- 
jection point is n s ~ 10 13 /cm 2 at V gs — 1.0 V if R s = 0.5 
ll.mm. 

The drain current and the transconductance are shown 
in Fig 2(b) as a function of the gate bias. The intrinsic 
transconductance peaks at g m ~ 1.7 S/mm. In addition, 
there is an inflection point at the gate voltage which cor- 
responds to n s = no at the source injection point. The 
inflection points are indicated by arrows. Experimental 
g m values as high as 0.8 — 0.9 S/mm have been reported 
recently. The phonon model thus predicts that as R s 
in GaN HEMTs is lowered, the g m will become more 
'peaky'. The inset (Fig 2(c)) shows saturated drain cur- 
rent as calculated from the phonon model (Eqs. 4, 3) 
in comparison with constant saturation velocity models. 
The phonon model follows a ^/n7 dependence at high 
densities (for example ~ 2 A/mm at n s ~ 2 x 10 13 /cm 2 ), 
and therefore explains experimentally observed currents 
in GaN transistors. 

The calculated current gain cutoff frequencies are 
shown in Fig 2(d) for L g = 50 nm using C g d = 2 
pF/cm. The gate dependence of follows that of g m , 
replete with a maximum near threshold voltage, inflec- 
tion points, and sharp drop at high gate voltages. Severe 
degradation of fr and 'smoothening' of the sharp fea- 
tures occurs with increasing R s and R^. To study the 
gate length dependence of /t, we choose the peak f T for 
each gate length and plot it against L g in the inset Fig 
2(e) for tb ar = 5 nm for the intrinsic case (fo = 1/2ttt c 
as n s — » 0), with C g d-, and with access resistances. A 
fr • L g = 23 GHz./im line corresponding to a constant 
electron saturation velocity of v sat ~ 1.4 x 10 T cm/s is 
plotted for comparison [7]. The phonon model predicts 
a slower increase in fa with scaling than the velocity 
saturation model for small gate lengths. 

The model presented here is able to reconcile the three 
peculiarities of GaN transistors and explains the current 
drive, g m , and fa behavior with gate bias. The limitation 
of the model as presented is that it assumes single sub- 
band occupation, which can be violated for high 2DEG 
densities. In such a case, the model can be extended by 
summing the current drive in each subband and incor- 
porating intersubband scattering. Effects such as double 
peak features in g m — V gs plots are then expected to ap- 
pear as higher subbands are populated. 
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